ROMA-1395/00 
SHEP 01/10 
hep-lat/010400^ 



K — > TTTT Decays in a Finite Volume 



C.-J.D. Lin% G. Martinelli^ C.T. Sachrajda'^ and M. Testa'' 

"'Dept. of Physics and Astronomy, Univ. of Southampton, 
Southampton, S017 IBJ, UK 

^Dip. di Fisica, Univ. di Roma "La Sapienza" and INFN, Sezione di Roma, 
Piazzale Aldo Mow 2, 1-00185 Rome, Italy 



Abstract 

We discuss finite-volume computations of two-body hadronic decays below the 
inelastic threshold (e.g. K — > tttt decays). The relation between finite- volume ma- 
trix elements and physical amplitudes, recently derived by Lellouch and Liischer, is 
extended to all elastic states under the inelastic threshold. We present a detailed com- 
parison of our approach with that of Lellouch and Liischer and discuss the possible 
limitations of the method which could arise due to the presence of inelastic thresh- 
olds. We also examine a standard alternative method which can be used to extract 
the real part of the decay amplitude from correlators of the form ( [TIttttT-Lv/K ] | ) . 
We show that in this case there are finite-volume corrections which vanish as inverse 
powers of the volume, which cannot be removed by a multiplicative factor. 
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1 Introduction 



The theory of non-leptonic weak decays represents a major challenge for particle physicists 
because of our inability to control the strong interaction effects. Lattice QCD offers a 
natural opportunity to compute the corresponding matrix elements from first principles. 
A number of theoretical questions however, needs to be answered before such computations 
can be performed with sufficient precision to be phenomenologically relevant. The main 
difficulties are related to the continuum limit of the regularized theory (the ultra-violet 
problem) and to the relation between matrix elements computed in a finite Euclidean space- 
time volume and the corresponding physical amplitudes (the infrared problem). The ultra- 
violet problem, which deals with the construction of finite matrix elements of renormalized 
operators constructed from the bare lattice ones, has been addressed in a series of papers 
0] and we will not consider it further in this work. The infrared problem arises from two 
sources: 

• the unavoidable continuation of the theory to Euclidean space-time and 

• the use of a finite volume in numerical simulations. 

An important step towards the solution of the infrared problem has recently been achieved 
by Lellouch and Liischer ||^ (denoted in the following by LL), who derived a relation 
between the K ^ titi matrix elements in a finite volume and the physical kaon-decay 
amplitudes. This relation is valid up to exponentially vanishing finite-volume effects. 

In this paper we present an alternative discussion of boundary effects and the LL-formula, 
based on a study of the properties of correlators of local s-wave operators Q We derive 
LL-formulae for all elastic states under the inelastic threshold, with exponential accuracy 
in the quantization volume. A consequence of our derivation is that the relation between 
finite-volume matrix elements and physical amplitudes, derived by Lellouch and Liischer 
for the lowest seven energy eigenstates can be extended to all elastic states under the 
inelastic threshold. Our analysis also demonstrates explicitly how finite volume correlators 
converge to the corresponding ones in infinite volume. 

Our approach is based on the property of correlators of local operators which can be 
expressed, with exponential accuracy, both as a sum or as an integral over intermediate 
states, when considering volumes larger than the interaction radius and Euclidean times 
< t ~ L. A possible source of concern is therefore the (practically relevant) situation in 
which the quantization volume allows only very few (perhaps two or three) elastic states 
under the inelastic threshold. In this case we will show that it is not possible, in general, 
to derive the LL-formula from the behaviour of correlators of local operators. At first sight 
this may appear to be a weakness of our approach compared to that of ref. [Q. We will 

^s-wave operators are those which, when appHed to the vacuum in a finite volume, excite only states 
with an s-wave component. They include the finite- volume restriction of all infinite- volume rotationally 
invariant local operators. 
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argue however, that the vahdity of the Liischer quantization condition for the energy levels 
of two-particle states in a finite volume p also requires the volume to be sufficiently large 
for the Fourier series to be equal to the infinite- volume energy integrals, up to exponential 
corrections. Since the derivation of the LL relation in ref. @] relies on this quantization 
formula, we conclude that the conditions on the volume for the applicability of this relation 
are equivalent in the two approaches. This issue is of considerable practical significance. 
For the foreseeable future, in actual lattice simulations of K it n decays there will be very 
few elastic states below the inelastic threshold, E^h (perhaps even only two such states). 
It is therefore necessary to examine the precision with which the integral over the energy 
is approximated by the sum over the elastic finite- volume energy levels below Eth. In this 
paper we discuss this question theoretically. 

In spite of the limitations discussed in the previous paragraph, we will argue that for an 
important class of dynamical situations the proximity of the kinematic inelastic threshold 
does not necessarily invalidate the LL-relation and that these situations are likely to include 
the two-pion system in physical kaon decays. A more quantitative answer will have to await 
results from detailed numerical simulations in the future. The number of states below Eth 
can be changed by varying either the masses of the mesons in a given volume (and hence 
considering unphysical decays) or the size of the volume. Thus one has the opportunity 
of investigating the importance of finite-volume corrections in exclusive two-body decays 
quantitatively. 

Another important consequence of our approach is the demonstration that it is possible 
to extract K ^ ttit amplitudes also when the kaon mass, mx, is not equal to the two- 
pion energy, i.e. when the inserted weak Hamiltonian operator carries a non-zero energy- 
momentum. Such amplitudes may be useful, for example, in the determination of the 
coefficients of the operators appearing in the chiral expansion. 

We also examine an alternative method often used to obtain information about K ^ nn 
decay amplitudes. This method is based on the evaluation of the correlation function 
( \7i^(ti)n_^(t2)'Hw{0)K{tK) I ), where Ti^it) denotes the three dimensional Fourier trans- 
form of the pion's interpolating field tt{x) f\ Taking the times ti^2 to be large and positive 
with ti ^ t2 and tx to be large and negative, we obtain the real part of the physical 
K 717T decay amplitude (see below), up to power corrections in the volume. These 
power corrections cannot be removed by multiplication by the LL-factor. 

The remainder of this paper is organised as follows. In sec. ^ we start by recalling some 
basic facts about infinite-volume Euclidean Green functions related to kaon decays. In 
particular, we stress the result derived in ref. |^ that in Euclidean space one obtains the 
average of the matrix elements into in and out two-pion states. We expand the discussion 
of the similarities and differences between the Euclidean- and Minkowski-space calculations 

^In the introduction we denote by K and tt interpolating fields with the quantum numbers of these 
particles. In subsequent sections we reserve K and tt to denote the corresponding particles and introduce 
a different notation for the interpolating fields. 
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in appendix A, where we demonstrate how the physical K ^ tttt amphtude (for which the 
two-pion state is an out state) is recovered in Minkowski space. 

We then proceed to discuss calculations in a finite volume. In section ^ we present a 
heuristic discussion concerning correlators in a three dimensional periodic cubic volume. 
In sec. ^ the finite volume formulae are derived and a detailed comparison of our approach 
with that of ref. is presented. In this section we also examine the conditions on the 
volume required for our derivation and in sec. ^ we analyze those needed for the Liischer 
quantization formula 0. 

In sec. |] we present an alternative method to obtain information about the kaon decay 
amplitude, based on the correlator ( \'7ig-{ti)7i-^{t2)'Hw{0)K{tK) \ ) introduced above. We 
show that in this case there are power corrections in the volume and explain why they are 
not removable by the LL correcting factor. We present our conclusions in section |^. 

There are three further appendices, introduced to clarify the discussion. In appendix B 
we demonstrate that it is indeed 3fJ{^} which is obtained from the {0\T[7niHwK] |0) 
correlator (3fJ{^} represents the real part of A) . In Appendix C we sketch the derivation of 
the Liischer quantization formula and some basic facts needed in the text. In appendix D 
we demonstrate that it is possible to determine the scattering phase S directly by computing 
four-pion correlation functions. 



2 Euclidean Green Functions and Physical Amplit- 
udes 

One of the main obstacles in the extraction of physical amplitudes from lattice simulations 
stems from the rescattering of final state particles in Euclidean space. The formalization 
of this problem, in the infinite-volume case, was considered in ref. |^ and is referred 
to as the Maiani-Testa no-go theorem. In this section we review the arguments leading 
to this theorem, elucidating the main differences between the Euclidean and Minkowski 
formulations of the theory. This allows us to fix the notation and to introduce several 
quantities which will be used in the following. The discussion of finite-volume effects is 
postponed to the following sections. 

We are ultimately interested in K ^ rni decays. Single particle states however, do not 
present any theoretical difficulty so, following ref. 0, in order to simplify the discussion 
while retaining the essential physical aspects, we eliminate the kaon and consider the Green 
function 

g{tu t2; q) = {0\ $ ^ (ti) $-,-(t2) J(0) I 0) , (1) 
where ti > ^2 > 0; J is a local operator which can create two-pion states with a definite 
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isospin, from the vacuum; 

%-{t)= I rf^xe-''?-^0(f,t), (2) 



and (f){x, t) is some appropriately chosen interpolating operator for the pion. For conve- 
nience the two-pion intermediate state is chosen to have total momentum equal to zero. 
From the continuation of the correlation function (|l]) to Minkowski space, using LSZ re- 
duction formulae one can determine the form factor 

out(vr(g)7r(-g-)| J(0)|0) , (3) 

which we assume to be the goal of some numerical simulation. In order to avoid unneces- 
sary kinematical complications, in the following we discuss the case of a Hermitian, local, 
scalar operator J, which excites zero angular-momentum states from the vacuum. The 
extension of the present discussion to correlators relevant for K ^ tttt decays is completely 
straightforward. 

In practice, in a numerical simulation the Euclidean correlation functions can be computed 
only approximately and so the continuation to Minkowski space is impossible. We there- 
fore have the problem of extracting the relevant physical information from the Euclidean 
correlation functions. 

The general expression of eq. (|lD in terms of matrix elements was derived in ref . 0] . For 

ti ^ ^2 > one has 

G{tut2;q) = ^[ out^(9)vr(-g) | J(0) |0) + MqM-q) I ^(0) |0) 



where 



+ 2E,Vg^{t2) I , (4) 
-P^ih) = ^E^(^"')l'^-^(0)l^)orout(r.|J(0)|0)e"(^"-^^')*% (5) 

n 

Eg = ^q^ + ml and = (0| 0(0, 0)|7r(g)) . (6) 

' Z and Eq can be obtained by computing the propagator of a single pion. We use the 
normalization ('7r(g')|7r(p)) = (27r)^2i?g(5^(g — p ). V represents the principal value of the 
integral implicit in the ^„ and the superscript "conn" implies that the corresponding 
matrix element is the connected one. Although the term in square brackets in eq. (^) is 
real and hence cannot be equal to the form factor in eq. (^, it is nevertheless a physical 
quantity which can be directly compared with experiment. To see this it is convenient to 
write 

3?(J) = ^[out(vr(g')7r(-g-')| J(0)|0) + in(7r(g)7r(-g)| J(0)|0) 
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,(7r(g)7r(-g)| J(0)|0)| cos5(2E,) , (7) 
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where S{2Eq) is the FSI phase shift for the two pions. Here and in the following we assume 
that the two-pion energy is below the inelastic threshold. Unfortunately, the extraction of 
3fJ( J) from the correlation function at large time distances is impossible in practice, since 
in this limit it is exponentially dominated by the last term in eq. (H). The conclusion of 
ref. was therefore that it is only possible to extract the form factor (7r(0)7r(0)| J|0), with 
the two pions both at rest. 

The continuation of eq. to Minkowski space is obtained by replacing t —>■ it in the time- 
dependent exponentials. At finite times therefore, the average in eq. also appears in the 
Minkowski space version of eq. (^). However, as shown in Appendix A, at asymptotically 
large time distances t2- 

2E,Vg-{h) I \ out(vr(g)7r(-g-) | J(0) |0) - in(7r(g)7r(-g-) | J(0) |0) 

SO that the terms in braces in eq. (1) combine to give the physical amplitude (including its 
imaginary part). In Euchdean space this is not the case. 



3 Finite- Volume vs Infinite- Volume Calculations: 
Three-Dimensional Case 

As pointed out in sec. 0, knowledge of 3?( J) together with the phase shift (5, is sufficient to 
reconstruct the physical amplitude. In infinite-volume the energy spectrum is continuous 
and the isolation and determination of J) turns out to be impossible. In a finite volume 
on the other hand, energy levels are discrete and the extraction of 3?(J) is possible in 
principle, provided that we are able to control the preasymptotic behaviour of the relevant 
correlation function at large time distances |^. 

We start this section by recalling the basic properties of the energy levels of two-particle 
elastic states in a cubic periodic box [^. In order to help clarify the arguments made later 
in this paper, a derivation of the relevant results is given in Appendix C. 

In a finite cubic volume, V = (where L is the length of each spatial direction), the 
allowed values of the "radial" relative momentum, fc, of a two-particle cubically invari- 
ant state with total momentum zero and with a non-zero s-wave component, satisfy the 
condition 

h{W, L)=n, (9) 
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where n is a non-negative integer Q and 

TT 

In eq. (jl^) 6{k) is the infinite-volume s-wave phase shift, q 
the center of mass energy W by 

W = + F . 

The function 0(g) is defined by 



(12) 
(13) 

In the following we will denote the energies (0) corresponding to each value of n hy En- 
Eq. (^ is derived under the simplifying assumption of an angular momentum cutoff: here 
and below we therefore ignore the contribution to eq. (^ from states of higher angular 
momenta. These could be included without difficulty in our discussion. Moreover we will 
consider 5 as a function of the total energy or of A;, whichever will be more convenient. 

The spectrum of cubically invariant two-pion states with total momentum zero is not 
exhausted by the solutions of eq. (P). There are states which correspond to g = |/| = |/'| 
such that I and /' G and are not related by a cubic rotation . It is shown in Appendix 
C that these additional states do not have an s-wave component, in contrast to those given 
by eq. d). 

Although the energy levels are those corresponding to a finite volume V , 6{k) in eq. (0) is 
the infinite-vohime phase shift. The derivation of this equation P], described in Appendix 
C, requires that V is larger than the two pion interaction region. In particular, in Appendix 
C it is shown that eq. (§) is equivalent to the condition that there is no distortion of the 
s-wave component of the two-pion wave function due to boundary conditions. We denote 
these wave functions by ^^^(x), where x is the relative position of the two pions. 

Our strategy is based on the analysis of correlators of local observables which admit two 
pions as intermediate states. For simplicity we take the two pions to have total momentum 
zero. We consider local observables which are rotationally invariant in infinite volume. As 
discussed in section HI this amounts to selecting the contribution of the s-wave component 
of the state \E'^J(x) to the matrix element (0| J(0)|-E'„). 

For n = 0, there are two solutions: one corresponding to fc = which is spurious, the other giving 
the Liischer relation between the finite volume energy and the scattering length . 



(10) 

kL/27T, and k is related to 



:iii 



tan 0(g) = ^ 



^oo(i;g2) 

where 
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In order to establish the relation between the finite and infinite-volume amplitudes, it is 
convenient to consider the two-point correlation function: 



C{t)= d-'x{0\J{x,t)J{0)\0), (14) 



V 



and its behaviour as the volume V becomes large. We have 

V 

|(0| J(0)|7r(pi)7r(p2))|' 



2(27r)3 
vr r dE 



J dEe-^'mmiTrn^E)]' J ^^6{p, + p,)6{E - E, - E,) (15) 
I ^e-^'mjmnn,E)\'k{E), 



where Ei^2 = \/Pi,2 + 



k{E) = \l—-ml. (16) 



On the other hand, in a finite, but asymptotically large, volume we have 

f d^x{0\J{x,t)J{0)\0)v = r^|(0|J(0)|7r7r,ri)y|V^"* (17) 

— > V [ dE pv{E)\{0\J{0)\7r7r,E)v\^e-^' . 



In eq. (P^), \n'K,n)v and \nn,E)v denote the finite-volume two-pion states at fixed n and 
fixed energy E, respectively. \7ni,n)v are classified according to the quantum number 
n defined in eq. (^. pv{E) denotes a function to be determined, which provides the 
correspondence between finite volume sums and infinite volume integrals. In many cases, 
for example in one dimension, pv{E) can be identified as the density of states at energy 
E. In section ||, after a discussion of the constraints imposed by the locality of J{x), we 
show that, also in the presence of cubic boundary conditions, pv{E) is given by, 

P^^^^ -dE- ^ ' ^^^^ 

with exponential precision in the volume. The expression in eq. ([T8|) is the one which one 
would heuristically derive from a naive interpretation of pv{E) as the density of states, cf. 
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eq. (^. There are however, some technical subtleties with such an interpretation which 
will be discussed later in this section. 



Comparing eqs. (|T5|) and (|T3), we establish the correspondence 



Ivrvr, E) ^ An ^^^^^y^ kvr, E)v • (19) 
Similarly, by considering correlators of operators which excite single particle states, we find 



\p = 0) ^ V2mV \p = 0)v . (20) 

From eqs. (|T^ and (^Up we readily obtain the amplitudes of the effective weak Hamiltonian 
in terms of the corresponding finite- volume matrix elements: 

\{nn,E\nwmK)\' = S2nW'P^^^^)^ |y(7rvr, E|7^^(0)|ir)v^P . (21) 

k{E) 

In order to obtain the physical amplitude we set E = tjik in the above equation and find 

3 



\{nn,E = mK\nwmK)f = SirV^ {q<j)'{q) + kS'{k)}k=k^ (^^^ |y(7r7r, E|7iiy(0)|K) 



|2 

v\ 

(22) 



where 



2 

m 



k^^^^-ml. (23) 

Eq. (p^) is the same as the formula derived by Lellouch and Liischer [Q. The additional 
factor of V"^ on the r.h.s. of eq. (^) is present because we have used the Hamiltonian 



density rather than the Hamiltonian. We stress that the relation (|21| ) is valid for the 
matrix elements of any local operator with any momentum transfer (below the inelastic 



threshold) and hence the matching condition tjik = E is not required. Although eq. ( p2D 
appears to be equivalent to the corresponding equation in ref. Q there is an important 
difference in the two derivations. The result of ref. was obtained at a fixed value of n 
and at a fixed volume V, tuned to have mx = En, with n < 8. Eq. (pi]), on the other hand, 
has been derived at fixed energy E for asymptotically large volumes V. This implies that 
as — cxD we must simultaneously allow n — *■ cxd. A question which arises naturally at 
this point is what is the relation between the two approaches? The answer requires a more 
detailed discussion, developed in the following section, where it will be shown that the 
constraints of locality allow us to establish eq. (|18D with exponential accuracy for elastic 
states under the inelastic threshold. 

We end this section with a discussion of the subtleties mentioned above, concerning a 
rigorous interpretation of pv{E) as a density of states. For purposes of illustration we 
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discuss the issues of completeness, normalization and degeneracy for the free theory. In 
this case the eigenf unctions are plane waves and we write the completeness relation directly 
in the subsector of zero total momentum 

piPn-{x-y) 

J2^^ = 6{x-y) x,yeV, (24) 

{Pn} 

where, as above, x denotes the relative position of the two particles. Since we will consider 
the insertion of this completeness relation between local, rotationally invariant operators, 
the only relevant component of the plane waves is their s-wave projection: 

VvW-4=/^.*-=4=^. (25) 

with Pn = |pn|- The component of the completeness relation which is equivalent to the 
identity when inserted between s-wave operators reads 

5^^n^p„(r)^p„(r')=X., (26) 

where Ig is the s-wave projection of the identity {5{x — y) appearing in eq. (p^) and Vn 
denotes the number of plane waves with the same p„, i.e. the same energy. The sum over 
states in eq. (p6|) is now a sum over the different energy values En rather than over the 
different momenta as was the case in eq. (plf). 



Eq. ( pq) shows that the factor Vn should not be considered as a degeneracy in the s-wave 
projection (it corresponds always to the same wave function); it should instead be viewed as 
a normalization factor. Thus the correctly normalized s-wave functions, when considering 
the completeness relation between rotationally invariant states, are 

t ( \ r I ( \ f^sinpnT 

fp^r) = V^n^pAr) = ^J ■ (2^) 

It is now straightforward to deduce the relation between the finite volume s-wave function, 
fp„{r) and the corresponding wave function in infinite volume: 

^ i^)(r) ^ ^ ! ^e^P- = ^ ^^^^^ (28) 
"^'^ V^(2^yc47r ^/(2^ pr - 

The relation is obtained by comparing eqs. (pTl) and (El): 
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Recall that Un in eq. ( p6D is the number of different ways in which a given n"^ can be 
realized using vectors with integer components. For large n, Vn is rapidly and irregularly 
varying with n and therefore the presence of a factor of Vn poses a difficulty to a rigorous 
interpretation of dn/dE as the density of states 0. However, what actually happens is that, 
whereas both factors pv{E) and |(0| J(0)|7r7r, i?)yp in eq. (p!7| ) are rapidly varying with 
energy, their product has a smooth behaviour. This can be seen by a computation of the 
correlator of a rotationally invariant observable, such as : (j)'^{x) : (where is a scalar field). 
For J(x) =: (f)'^{x) : eq. ([T7|) becomes 



(30) 



d^x{0 \ : 0^(x) : : 0^(0) : \ 0)v 



V 

f c^^^a; V |(0| : 0^(0) : |pi,P2)y|' e'^^'+P'>^ e-^''-^^''-^^ 
^ Jv ^ ^ 

Pl,P2 



2 
V 



E 



(01 : ,#.^(0) : |p, -p)v 



V$:i(0| : 0^(0) : |nMV-- = ^5:4 

oo 

■ /" dE pv{E) |(O|:02(O):|n)vpe-^* 



z/„ e 



In eq. (0) we have used the free-field result 

(O|:02(O):|p,-p-)v^ 



VEr, 



so that 



|(O|:0^(O):|n)p 



\/2 ''^ 



Therefore the correspondence in eq. (|T^) is given in the free theory by 



{0\:m--\PuP2)\' = \{0\:m--\n)\'El '^'^^^ 



" 27r2n2 



(31) 



(32) 



(33) 



In deriving eq. (|33D we have used the fact that in the free case the spectrum of states, and 
hence the solution of eq. (||), is given by 



2 / Lk\ 2 
^ =t,2^J =^ 



(34) 



*Wc acknowledge many discussions on this point with L. LcUouch and M. Liischer. 



10 



It is easy to demonstrate that 



<f>'{n) = , (35) 

and hence that 

|(O|:02(o):|p,,p2)|' = 4, (36) 

which is the correct result for the free-field theory with covariant normalization in infinite 
volume. 

Following this discussion, the last line of eq. (|5D|) is equal to 

oo 

vldEpy{E)\{0\:<p\0):\n)y\'e-^' = ^^l^^dE^e-^' (37) 



in perfect agreement with the infinite volume expression. 

This example shows that, in the free case, the meaning of the LL factor is very simple. It 
removes the degeneracy factor (i/„) due to the fact that we extract from the correlator the 
contribution of the n^^ energy level and restores the appropriate factors to normalize the 
states in the infinite volume limit correctly. 



4 How large is a large volume? 

In this section we study finite-volume effects in a toy mathematical example as well as 
in non-relativistic quantum mechanics and relativistic field theory. The common feature 
and key ingredient in all these cases is the locality of the operators whose matrix elements 
are being computed. Within our approach we will be able to show that, in the quantum 
mechanical case, the LL formula is true for all levels containing an s-wave component, up 
to exponentially small corrections in the volume. This conclusion is compatible with a 
stronger unpublished result by Liischer |§] , who finds that the LL relation is indeed exactly 
valid for the same levels. The discussion given in ref. [Q, however, is specific to non- 
relativistic quantum mechanics and cannot be readily generalized to relativistic quantum 
field theory. 

In relativistic field theory the situation is more complicated for two basic reasons. Firstly, 
the interaction range is not precisely defined, so that all results are valid at most with 
exponential precision. The second and more substantial difficulty comes from the presence 
of inelastic production, which, as shown in sec. ^ below, could also affect the validity of 
the quantization condition in eq. (||). 
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For the discussion below, it will be useful to recall the general relation between a discrete 
sum and the corresponding integral, which follows from the Poisson formula |0 



oo 



Y,S{x-n)= ^ e^2— . (38) 

n=— oo n=— oo 

Multiplying both sides of eq. ( ^8|) by a function f{x) and integrating over x, we get: 

poo 

E /(^) = E / (39) 

n=—oo n=~oo 

/oo poo 
f(x)dx + J2 / /(a:)e2™^cia;. 
-oo 

Eq. (^) shows that an integral is equal to a discrete sum up to correction terms, given 
by the sum over terms with n ^ 0. In the following we will show how to control these 
corrections in order to make contact with the LL formula. 



4.1 A Toy Example 

Consider the simple one-dimensional example of a C°°-function, f{x), vanishing with all 
its derivatives outside a compact support which is entirely contained inside the "finite 
volume" —L/2 < x < L/2. This is the statement of locality in our toy model. f{x) can be 
expanded as a Fourier series: 



L 

n 



where 



and 



Pn = j^n (41) 



r.L/2 

f{p)= I dxf{x)e~'P''=l dxf{x)e-'P\ (42) 



' -L/2 

On the other hand, f{x) can also be written as a Fourier transform 



f{x) = ^l dpf{p)e'P\ (43) 
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Whereas eq. (|i3| ) reproduces /(x) exactly, eq. (|iO| ) coincides with /(x) only in the interval 
—L/2 < X < L/2] outside this interval f{x) is replicated periodically on the x-axis. If we 
now consider a second function g{x) with the same properties as /(x), we have: 

dxg*{x)f{x) = Y,~g*{p,)f{p^) = / dp~g*{p)hp) . (44) 

Eq. (^if) is an example of the exact equality between an integral and a discrete sum without 
any correction terms. It shows that, for local observables, the "finite volume correlator" 
T J2n9*(Pn)fiPn) IS Gxactly equal to the "infinite volume correlator" 2^ f^°° dpg*{p)f{p), 
as long as the support of the observables is entirely contained inside the "quantization 
volume" L. In the next section we will show how to derive similar relations to the LL- 
formula from equations such as (H). 



4.2 Finite- Volume Effects in Non-Relativistic Quantum Mechan- 
ics 

We now extend the previous considerations to the case of a quantum system, starting 
with the non-relativistic quantum mechanical case. In preparation for the relativistic 
case presented below, we use the formalism of second quantization also in the quantum 
mechanical case and introduce interpolating operators 0^(x, t) and 0(x, t) that respectively 
create and annihilate the physical particles. Consider two "local" , zero angular-momentum 
two-particle operators TZ{X,t) and S{X,t) defined by 

n{X,t) = ^ j d'^xR{r)(P^ (^X + ^,t^ (f)^ (^X , where r = |x| , (45) 

with a similar definition for S{X,t). The functions R{r) and S{r) defined in eq. (^5[) have 
support contained inside V, thus ensuring locality. Under this hypothesis, we now show 
that the matrix elements of TZ and S in the finite volume V are proportional to those in 
the infinite volume, with a constant of proportionality given precisely by the LL-factor. 

Consider the zero momentum two-point function 

gv{t)= [ d^x{o\s\x,t)n{o,o)\o) . (46) 

Jv 

At t = we have 

GviO) = J d'xSir)R{r) = j d^xS{r)R{r) = g^{0) , (47) 

where the second equality comes from the hypothesis that S{r) and R{r) have support only 
in V. Eq. (^) shows that we can analyze Gvi^) both in terms of finite volume and infinite 



13 



volume intermediate states. Inserting a set of complete states in the infinite volume we 
have 

POO 

6?oo(0)= / dEs*{E)r{E) (48) 

J E 

where E is the lowest energy in the vrvr channel, £' = or 2m7r in the nonrelativistic and 
relativistic cases respectively and, for example, 

s*{E) = j d^xS{r)m%-'^''"'\r) . (49) 

In eq. ( ^9|) we have denoted by \l'g ^"■'"^'> (j.'^ ^j^g center of mass, zero-angular momentum 
wave function of two pions with energy E. 

We now turn to the correlation function in a finite volume, and insert a complete set of 
energy eigenstates: 

^^(0) = y ^( |5t(0, 0) |n) (r^|7^(0, 0) | ) , (50) 

where \n) denote the energy eigenstates with zero total momentum and with the normal- 
ization {n\n') = 6n,n'- Since the local operators under consideration depend only on r, the 
sum in eq. ( pOD receives contributions only from the zero angular momentum component 
of the intermediate states. In other words we have, for example 

(0|5^(0,0)|r2) = [ d^xS{r)'^^p{x) oc s*{En) . (51) 



V 



As shown in Appendix C, the states |n) are precisely those described by eq. (|^). We may 
therefore rewrite eq. ( ^OD as 

Gv{0) = vY,s*{E^y{En)cv{En), (52) 

n 

with a set of non-negative coefficients cv{En), in terms of which we have 

\s{En)\' = ^^A{0\S\6Mn)\' . (53) 

Cv{En) 

It is our aim to evaluate the coefficients cv{En) and to demonstrate that eq. ([53|) is precisely 
the LL-relation. 

The previous discussion shows that at t = the correlation function Qv{^) defined in 
eq. (|46|) coincides with the corresponding one in the infinite- volume limit, without any 
corrections. We now investigate what happens at a generic value of t. In this case we get 

gy(t) =V ^{0\S^{0,0)\n){n\'JZ{0,0)\0)e-^"' = j d^xS{r,t)R{r) (54) 
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where 

-m 



5(r,t) = e"^'5(r), (55) 



with a similar result for i?(r, t). In eq. (|55D H denotes the first quantization, center of mass 
Hamiltonian (differential) operator acting on two-particle wave functions, whose eigenval- 
ues are precisely the {En}- From eq. (|55| ) it follows that S{r, 0) = S{r) and i?(r, 0) = R{r). 

The key point for the identity of finite and infinite volume correlators at t = 0, as seen 
from eq. (^), is the locality of S{r) and R{r). On the other hand, for a generic t, eq. (^Sj) 
shows that ^(r, t) is obtained from S{r) through a diffusion process. It is well known that, 
even if S{r) has a compact support, S{r,t) will be non-zero over the entire infinite volume 
for t ^ 0. In general, therefore 

Gvit) ^ Qooit) . (56) 
However, diffusion is a "slow" process, and for a range of t ~ L we still have 

6^y(t) = 6^oo(t) + 0(e-"-^) , (57) 



while for t > m^L^ we expect Qvit) to be very different from ^oo(^)- Eq. (|57D implies that, 
with exponential accuracy in the volume, 

POO 

y Vs*(K)r(K)cy(K)e-^"*= / dEs*iE)r{E)e-''r (58) 

We can now use eq. (^) to transform the r.h.s. of eq. ( p8|) back into a sum again. Eq. ( ^8]) 
implies that 

oo oo 
J2 HHE, L)-n)= ^^^nh(E,L) ^ ^gg^ 



where h{E,L) is defined in eq. (|10D . 
Combining eqs. (|58D and (p^) we obtain 



r rfE.*(E)r(E)e-^* = ^ - Q(i^, t) , (60) 

n \ dE lEr, 

where the observable S and TZ have been chosen to vanish sufficiently rapidly in E to 
avoid problems with the spurious solution with k = (see below). In eq. (^), we have 
introduced 



POD 

Q{L,t) = Y, ciEs*(E)r(E)e-^*e*2-"^(^'^) 
1^0 -^^ 



(61) 
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and used the fact that the energies are non-negative. 

A general property of any quantization condition is that for fixed E 

h{E,L) — > oo . (62) 

In our case, with cubic boundary conditions, we have 

h{E,L) — > Lk{E). (63) 

L— >oo 

Eq. (^) is a statement that as the volume is increased, the value of n corresponding to 
any fixed energy also increases. 

For n = 0, a possible complication may arise due to the spurious solution with = 0. 
This can be easily avoided however, since it is possible to choose local operators for which 
the threshold behaviour of s*{E) r{E) can be made to vanish arbitrarily quickly. This 
can be achieved by taking appropriate combinations of local observables and their time 
derivatives and recalling that any time derivative multiplies the matrix element by the 
energy according to 

{0|5^(0,0)|n) = i{0\[H,S\0,0)]\n) = -iE{0\S\6,0)\n) . (64) 



Using the techniques of asymptotic analysis |jTO|, and taking into account eq. (|6^) , we 
deduce that the behaviour of the correction term Q{L,t) at large L is dominated by the 
critical points of h{E, L), i.e. the points at which h{E, L) has a vanishing derivative or the 
points where h{E,L) and s*{E) r{E) are not continuous or differentiable. In the present 
case the only critical point of h{E, L) occurs at threshold, i.e. for k ^ 0. This means 
that the behaviour of Q{L,t) at large L is power-like, with the power determined by the 
threshold behaviour of s*{E)r{E) e~^^ . These power corrections depend on the time as 
e~^^ ~ g-2m,^i (^multiplied by powers of t). They therefore only affect matrix elements 
computed at energies close to threshold. Indeed, as discussed above, by an appropriate 
choice of local observables it is possible to make the product s*{E)r{E) at threshold ar- 
bitrarily small and consequently to make the leading power in the 1/L expansion in Q 
arbitrarily high. Thus 

Vj2s*iE^)r{Er.)cviE^.)e-^-^ = '^^^^^(if - Q{L,t) , (65) 

and since we are interested in amplitudes at fixed energy E (between E and the inelastic 
threshold), in the limit L oo, n oo, we have, with exponential accuracy, 

CV-l^n) = ^^ 9h(l,L) \ (66) 
^ I dE \e„ 
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for all states with energies above E. Eqs. (|53| ) and (|66[), together with the kinematic factors 
appropriate to the chosen normalization, establish the LL-type relations of the kind given 
in eq. (0), with exponential accuracy. 

We conclude this subsection with some brief remarks about the matrix element with the 
two-pion state at threshold. A priori, the threshold matrix element can be affected by power 
corrections, although our derivation does not necessarily imply that this is the case, since we 
are unable to control the sum over I in eq. (|6ll). Threshold states are generally problematic 
and require a separate discussion. For example, in the case of an attractive potential, 
their energy in a finite volume is lower than E. We are not able therefore, to derive 
a relation between the finite-volume matrix element at threshold and the corresponding 
infinite-volume amplitude which we can demonstrate is free of corrections which vanish as 
inverse powers of the volume. 



4.3 Finite- Volume Effects in Relativistic Quantum Field Theory 

In relativistic quantum field theory (RQFT) the previous considerations must be modified 
since strict localization of the support is not possible and therefore the state 

|J)^J(0)|0) (67) 

is only approximately localized around the origin. The probability of finding particles at a 
distance r away from the origin never vanishes, although, as a consequence of clustering, it 
decreases exponentially like e~"^^^ . There is a further point that needs to be considered. 
The correlator appearing in eq. ( P7| ) is evaluated at equal times. In field theory such an 
equal-time correlator is generally UV divergent, due to the short distance singularities 
in the product of two local operators. A natural regulator for such divergences is the 
(Euclidean) time. The discussion of the quantum mechanical case can then be generalized 
to field theory. We therefore conclude that the equality 

Gvit)=gooit) (68) 

is also valid in RQFT with exponential accuracy. 

For correlators of local operators, the most important difference from the non-relativistic 
quantum mechanical case is due to the existence of the inelastic threshold, Eth- This 
implies that in any fixed finite volume, there are only a limited number of elastic states 
with energies below Eth- Therefore the truncated sum over the En (truncated at Eth) 
reproduces the corresponding integral (cut-off at the same energy) only up to corrections 
which vanish as L — »• cxd. However, the number of such states increases with the volume 
and the arguments presented in section |4.2| can readily be generalized to the RQFT. In 
particular we conclude that eq. ( pB| ) is still valid with exponential accuracy. Our derivation 
of the LL formula demonstrates that, in principle, it is possible to improve the precision 
by working at fixed energy on increasing volumes (and hence increasing n). 
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The discussion above does not imply, however, that the corrections are small for any vol- 
ume. In order for these corrections to be really negligible, we must work in volumes which 
are sufficiently large that the sum over energies is well approximated by the corresponding 
integral. It is of course impossible to give an estimate of the error on the LL formula in 
lattice simulations when the volume is such that there are only a few elastic states (2 -^ 3) 
below Eth- We will argue in the next section that in such cases there are related difficulties 
in the derivation of the quantization condition (^. 



5 Quantization of Two-Particle Energy Levels in 
Quantum Field Theory 

In this section, we re-examine the quantization condition given in eq. @ within the frame- 
work of Relativistic Quantum Field Theory (RQFT). Although eq. (j^ has been discussed 
in detail in the quantum mechanical case, and its extension to RQFT can also be found 
in ref. 0, to our knowledge the explicit, non-perturbative treatment given below is new. 
This discussion helps to clarify what happens when only few elastic levels exist below Eth- 

The concept of a wave function in RQFT is an approximate one. In this case, the object 
closest to a wave function is the Bethe-Salpeter (BS) wave function. We discuss the infinite- 
volume case first. For an incoming state with total momentum zero, the BS wave function 
is defined as 

%^{x) = ( I r[ 0(f , t) 0(0) ] I n{p) 7,{-p) )i, . (69) 

For simplicity we present our discussion for t = [|. We now derive the relation between 
$p'(x) = $p-(x, 0) and the scattering phase shifts in the infinite-volume limit. We introduce 
a complete set of intermediate states 

%<x) = 5^(O|0(f,O)|n)(n|0(O)|7r(p)7r(-p))i„ (70) 

n 

= i^q) \m\Ap)A-p)hne"^-+I{x) , (71) 

where X(x) represents the contribution coming from inelastic intermediate states with more 



than one particle. We separate the connected and disconnected contributions to eq. ( ffl]) 
in the standard way: 



%^{x) = Z e^^- + / ^ ( ^i^) I m I ^iP) <-P) )'^^'' + m . (72) 



^There are no ultraviolet problems, even if (j){x) is a composite field (as is the case in QCD) because 
we always choose x 7^ 0. 

^In order to simplify the presentation we will treat the particles as distinguishable. 
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Using 



imo) I MP) A-P) >- - --i^^, - f ^^3^ . (^3) 

we obtain 

Ml) - Z (e-- + / S JEjE^^W;^r) + ^(--) ■ (^4) 

For large values of x the inelastic contribution can be safely neglected. Moreover, in 
eq. (0), the integral over g is dominated by the on-shell contribution [E^ ^ E'p), for which 
becomes the physical on-shell amplitude M.[j) q). Indeed M. has no singularities 
for 2Ep below the inelastic threshold. To simplify the discussion, we consider a situation 
in which all the phase shifts except the s-wave one (denoted by 5{p)) are negligible. Under 
this hypothesis we have 

M{p-^ q) = — ^ (e'''^^^ - 1) . (75) 
I p 

From eq. ([74D, we then get 



(1 pipr \ 

2i ^ pr J 

Using the asymptotic behaviour of the s-wave projection of e*^'^ 

sinpr 



\s—wave 

pr 



+ . . . , (77) 



we finally arrive at 



$p(f)|™e = -Lmpr + ^.{e^''^'''> -l)e'A 
pr \ 2i ^ J 

= sin {S{p) + pr) + . . . . (78) 

pr 

We are now in a position to discuss the quantization of energy levels in a finite volume. 
For our purposes it is sufficient to consider the simple case of a spherical finite volume of 
radius R with vanishing fields on the boundary []. For sufficiently large R, the asymptotic 



expression in eq. (|78D to hold, we obtain the quantization condition 

nn - 5{p) = pR , (79) 



^We return to consider the quantization condition in a cubic volume in some detail in Appendix C 
below. 
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which is the analog of eq. (which was derived for a cube) on a spherical volume. 



The central question for our discussion is whether the quantization condition in eq. ( [79D 
holds when the volume is so small that only a few levels exist below Et^. In such a situation 



we show that in general it is not possible to follow the steps that lead from eq. (^) to ( [79D 
from which the quantization condition is obtained. 

As previously discussed, the validity of the Liischer quantization condition requires (the s- 
wave projection of) $p (x) to be undeformed by the presence of the boundary. This implies 
that, up to terms which vanish exponentially with the volume, it is the same function of 
X (for xeV) when expressed as a sum over either the finite or the infinite-volume energy 
eigenstates 

= 5^(0|.#.(0)|n>(r#(0)|7r(p-)7r(-p)>i.e■«•■•"^ (80) 

n 

where are the 3-momenta appropriate to the given boundary conditions. 

Both the 5-function and the principal part of the pole at Eg = Ep in eq. (^) are essential 
to obtain eq. ([76|); they separately contribute to both the outgoing and incoming waves and 
it is only in their sum that the principal part cancels the incoming wave and doubles the 
outgoing one. It follows that the scattered wave is purely outgoing, as expected. We remind 
the reader that for this derivation of the quantization condition to be valid it is necessary 
to work at distances such that the integral is dominated by energies Eg ~ Ep, so that the 
variation of Ad with the energy can be neglected. Thus a necessary condition for the validity 
of the quantization formula in a finite volume is that the same situation is reproduced by 
the discrete sum in eq. (|^). Strictly speaking, in a fixed volume, it is not possible to 
separate the connected and disconnected terms and to parametrize the amplitude as in 
eq. (|73D , because no energy poles can be present in this case. These contributions and, in 
particular, the principal part of the integral in (|7^ must be well approximated by the sum 
over those elastic states for which we may neglect the variation of A4. Therefore, for the 
validity of eq. (^), the sum over n in eq. ( ^OD must be dominated by its elastic part. On 
the other hand, we expect that, even restricting the sum to one-particle intermediate states 
{\n) = \qn)), the matrix element {qn\(f>{^)\ ^(p) ^(~P) )m receives important contributions 
from the inelastic states when En > Eth.. This happens because there is no reason for 
the off-shell amplitude (n|(/)(0)| 7i{p) 7t(— p) )in to be dominated by the one particle term. 
We conclude, therefore, that in the case where only a few levels have energies below Efh, 
the sum over the finite- volume eigenstates cannot, in general, be dominated by the elastic 
contribution, and we should expect corrections to eq. ([791). 



An equivalent way of arriving at the same conclusion is to write the Schrodinger-like 
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eigenvalue equation which is satisfied by $p'(x), namely 

{ml - - El) %^{x) = I Md^-' + K{x) ^ /CeKf ) + /C(f ) , (81) 

where K,{x) = {ml — — E^) X(x). In the non-relativistic case the r.h.s. of eq. (|8T|) is 
simply V{x)^^{x) and is therefore well localized, as discussed in the previous section. In 
the field theoretical case, /Cez(^) + }C{x) is an exponentially localized function in x, since it 
is the Fourier transform of a regular function. The differential operator {ml — — El) 
applied to ^p{x) filters out the long distance contribution from the disconnected and pole 
terms, leaving only the short range contribution. Thus, for the quantization condition to 
be valid }Cei{x) + /C(x) must be completely localized inside the quantization volume. Once 
again this is true if and only if lCei{x) + lC{x) is the same function, when expressed as a 
Fourier integral or a Fourier series. However, while in the infinite-volume case there is a 
substantial contribution to ICei{x) + IC{x) coming from elastic states with energies below 
Eth, in a finite volume this contribution is limited to a sum over a few (perhaps even as few 
as two or three) elastic states and the Fourier series and integral cannot be expected to be 
identical. It is true that, in the sum over intermediate states, we still have an infinite sum 
due to the contribution of all the states above Efh, however, in cases in which inelasticity 
is important, this contribution is qualitatively different from the ones below Et^. 

We end this section by noting that there is a favourable dynamical situation in which 
the quantization conditions (^ (and (^)) remain valid and the LL correction formula 
is applicable even on volumes with only a few states below Eth- This is the case for 
weakly interacting particles, for which inelastic contributions are negligible at the energy of 
interest. Such a dynamical situation is effectively the same as in quantum mechanics, where 
it is sufficient for the interaction region to be only slightly smaller than the quantization 
volume. From phenomenological studies it is likely that such a situation describes, to a 
good approximation, the low energy dynamics of pions [|1^. 

6 Evaluation of the amplitudes: a different method 

In this section we analyze a different way to obtain weak K amplitudes in finite- 

volume simulations, one that is more closely related to the discussion in ref. and sec- 
tion 0. The method, which is currently used in numerical simulations, requires the eval- 
uation of four-point correlation functions of the form {0\7i^(ti)n_^(t2)Hw{0)K(tK) \ 0), 
where the times ti^2 are large and positive with ti ^ ^2 and tx is large and negative. We 
will show that from this correlation function one obtains directly the real part of the decay 
amplitude, up to corrections which vanish as inverse powers of the volume and which can- 
not be removed by a multiplicative correction factor. The connection of this calculation 
with the LL approach will be explained. 
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Since single-particle states do not present any theoretical complications, we proceed as 
in section ^, eliminating the kaon and considering the correlation function Q defined in 
eq. (|l]).When ti oo, eq. can be conveniently written as 



S{ti,t2;q) ^ T^^e-^"*^ 6^3(^2, g"), (82) 
where 

6^3(t,g1 = (7r(g)|<f-,-(t)J(0)|0), (83) 
for t > 0. From eqs. (§)-(0) we obtain 



G3{t,q) = ^e-^''' |out(7r(g)7r(-g)| J(0)|0)| cos5(2E,) +■■■ , (84) 



^1 



where the ellipses represent the principal value term defined in eq. (^). We remind the 
reader that the above equations are valid in infinite volume. 

If we could isolate the term proportional to exp(— i^^t) in the correlation function ^3, we 
would obtain |out( 7r(g )7r(— g ) | J(0) | 0)| cos5{2Eq). This can be achieved by working in a 
finite volume in the following way. Considerations similar to those in section ^ give 

00 

6;3(t,g) = VWE^v'(i^^n)e-(^"-^')* / ciEAoo(i?)e-(^-^«)*, (85) 

where the finite- volume spectral function Ay is defined as 

Av{En) = y(7r(g) I $_,-.|7rvr, n)y y ( vrvr, n | J(0) | ) . (86) 



Eq. ( P^D allows us to write 

00 

j dE A^{E)e~^''-^''^' = 

e-^^* |out(vr(g)7r(-g)| J(0)|0)| cos5(2E,)+ J dE A^''iE)e-^^-^^^\ (87) 



2m„ 

'z 



2E, 

2 



where P.V. denotes the contribution coming from the principal value of the integral. Using 
the Poisson summation formula in eq. (^), we obtain 



^3(t, q) e-(^--^')* I out( 7r(g)7r(-g ) | J(0) | ) I cos5(2E,) 



V^oo 2Eq 



+ J:^^'-'""-"'''+0(t) , (88) 
^ Pv{En) \LJ 
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where En is the closest approximation to 2Eg among the solutions of eq. (|]). The replace- 
ment of 2Eq by Efi, at fixed t, is valid up to power corrections in the volume. Due to 
the presence of pv in the denominator, each term in the sum in eq. (^) vanishes at large 
volumes (as inverse powers of the volume). Thus, from a numerical evaluation of Q we 
obtain 

AyiEn) = ^ I out(7r(g)7r(-g) | J(0) |0)| cos5{2E,) , (89) 

up to power corrections in the volume. Eq. (|89|) was suggested in ref. 0, under a smooth- 
ness assumption on the amplitudes. 

We end this section by explaining that the conclusion of the preceeding paragraphs is 
consistent with the results of ref. 0] and the discussion in the previous sections. Consider 
eq. (^). From the final-state interaction theorem, we deduce that the phases of the two 
factors are equal and opposite so that the right-hand side is equal to the product of the 
absolute values. By an additional measurement of the four-pion correlation function (see, 
for example. Appendix D), we can obtain |y(7r(g) | ^-g- Itttt, n)v/|, and hence determine 
|y( VTTT, n I J(0) I ) |. Correcting this finite- volume matrix element by the LL-factor we 
determine the infinite-volume matrix element up to exponentially small corrections. This 
is not the procedure which was followed in arriving at the real part of the matrix element 
in eq. (^) however, and this explains the different precision of the two cases. 



7 Conclusions 

In this paper we have discussed the conditions under which the amplitudes of two-body 
decays below the inelastic threshold, such as K ^ -k-k decays, can be computed in a 
finite volume. We have shown that the correction factor recently found by Lellouch and 
Liischer 0], which relates the finite- volume matrix elements and the physical (infinite- 
volume) decay amplitudes at zero momentum transfer for finite volume states up to n = 7, 
can be extended, with exponential precision in the volume, to all elastic states under the 
inelastic threshold and to any momentum transfer. 

We have examined the possible influence of inelastic thresholds both on these correction fac- 
tors and on the quantization condition for two-particle states (below the inelastic threshold) 
in a finite volume. Our conclusion is that, in general, the presence of inelastic thresholds 
requires us to work on sufficiently large volumes so that the sum over elastic states under 
the inelastic threshold can be approximated reliably by the corresponding integral over the 
energy. In lattice simulations in the foreseeable future the fulfillment of this condition will 
represent a formidable challenge. It is therefore reassuring that in hadronic sectors, such 
as the two-pion system m K -k-k decays, where the interaction is rather weak and in- 
elasticity develops fully only at energies significantly higher than the kinematic thresholds. 
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the finite-volume approach may be applicable starting from relatively small volumes. The 
quantitative implications for lattice computations of — > vrvr decay amplitudes will need 
to be investigated in numerical simulations. 

We have also examined a different and frequently used method to obtain the real part of 
the K ^ TTTT amplitude based on the correlator (0 | T[7nrT-CwK] \ 0). We find that in this 
case there are finite-volume corrections, which vanish only as inverse powers of the volume, 
and we have explained the connection with the approach of ref. @]. 

Nonleptonic weak decays will continue to play a central role in particle physics phenomenol- 
ogy in the coming years. We trust that this paper is a contribution to the theoretical 
framework which will underpin lattice simulations of K —>■ utt decays. Further work is 
needed in order to understand how to treat contributions from states above the inelastic 
threshold, particularly in relation to nonleptonic two-body i?-decays. For this important 
class of processes a huge amount of experimental data is becoming available and yet its 
interpretation is currently limited by our inability to quantify the strong interaction effects. 
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Appendix A 

In this appendix we show that Q(t,q), defined in eq. (|l|), is equal to the expression given 
in eq. also in Minkowski space. We will see that although the right-hand side of eq. (^) 
contains the average of the matrix elements into in and out two-pion states, the term 
containing the principal value integral is precisely that required to recover the full physical 
amplitude in Minkowski space, but not in Euclidean space. The aim of this discussion is to 
clarify the the similarities and differences between the Euclidean and Minkowski correlation 
functions. 

We start by inserting a complete set of out-states {\m)} in Q^, defined in eq. (|83|), and 
^ As before we restrict our discussion to two-pion intermediate states. 
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obtain 

Qsit, q) = Y,{27rf 6^'\pm) e-(^--^^)* {7r{q) \ $_,-(0) |m)(m | J(0) | 0) , (90) 

m 

where Pm is the momentum of the two-pion state labeled by m. 
The disconnected contribution to Gs is readily obtained using 

{n{q) I $_,-(0) \n{k)n{-k))t:^ = 2E, {27:f 6^'\q-k) y/Z , (91) 

so that 



^3^^^(^, Q) = ^ e-^-'* om{AQM-Q) I ^(0) |0) . (92) 

To evaluate the connected contribution we use the reduction formula 

out("^ 1 7r(gi)7r(g2) )in = 2E2 ont{m-iT{q2) 1 7r(gi) )in 

+-^(27r)V^)(p^-?i-g2)(-p' + m^) out(m |</)(0) | gi ) , (93) 

where temporarily we introduce general momenta qi and g2- {fn — 7r(g2)| represents the 
state containing the particles in state (m| with n{q2) removed. This contribution is absent 
if (m| does not contain such a pion. The scattering amplitude M.{qi,q2;m) is defined from 
the hmit as p = {Em-Ei, q2 ) goes on-shell, — > m^, 

(7r(gi) 10(0) I m)^°r = , ^\ = - V^t ^ ■ (94) 

/out _p2 + ^2+^^ [p^_E2-ie){n + E2 + ie) ^ ' 

The physical scattering amplitude for 7r(gi) + 7r(g2) ^ m is therefore given by 

(27r)^5(^)(p„ - ?i - 52) («Al(gi, ?2; m)) = out("^ |7r(5i)7r(g2) )in - out("^ \'J^{q\)'K{q-2) )out , 

(95) 

where we have used the fact that single-particle in and out states are the same. 
We now use 

^ i7r(5(po-^2) +7^^— , (96) 



7-1 • " V-i^U ' ' 7-1 ' 

Po-E2-te P0-E2 
where V stands for "principal value" , to write 



(7r(g0 10(0) I m)— --i^^M*- ^r,^* V-^ . (97) 
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The first term on the r.h.s. of ( p7[ ) gives a contribution to Q3 of 

MQ1MQ2) |J(0) |0)- out(vr(gi)7r(g2) |J(0) |0)] , (98) 
which, when combined with the disconnected term gives 

,(7r(gi)7r(g2)|J(0)|0) + out(vr(gi)7r(g2) | J(0) |0) 



^conn,A ^ VZ iE^t^ ^ 1 

^ 2E2 2 



^3 "2^;^ 2 



(99) 



Thus we see that even in Minkowski space we get the average of the matrix elements into 
in and owt-states. The remaining contribution to Q-^ comes from the principal value term 
in (13). 



In Minkowski space however, writing the result in the form of (^) is not very transparent. 
To recover the standard result recall that ^3 contains a factor of exp(— iii^m^) and that 

/ I^W^ /aw \conn \fZ M* \fZ M* ,^ 

l*-,-(0) = -^f—r— = - E„XE^-2E,-i^ ■ 



The sum over the intermediate states can be written as an integral over Efyi' As t — > 00 
this integral, which is sufficiently convergent, can be extended to the range (—00, 00) and 
evaluated by contour integration. For positive t we can close the contour in the lower 
half-plane and, since there are no singularities there, we obtain zero for the connected 
contribution. In this way we obtain the standard result 

G'''''"ih,t2;q) = -|-e-^^'(*^+*^) out{niq)ni-q) \ J(0) |0) . (101) 



In Euclidean space it is not possible to treat the connected part in this way. Instead of the 
factor of exp{—iEmt), Q3 now has one of exp(— i^^t) which, at large times t is dominated 
by Em — 2m^ i.e. the lowest energy state consistent with three-momentum conservation. 
It is therefore not possible to perform the energy integration by contours. 



Appendix B 



In this appendix we present a derivation of eq. (184| ) based on elastic unitarity. The starting 
point once again is the quantity 

g,{t, q) = {n{q) |<l>_,-(t) J(0) | 0) , (102) 
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defined in eq. ( |83| ) of section ^. From tlie discussion in sec. Q and appendix A we know 
that ^3 is given by: 

r?. ^ r . 1 

+ 



^3(t,g) = ^e-^^\,,(7r(g)vr(-g)| J(0)|0)x [l-^A<*(2i?,)Ar(2i?J 



V I ^f^{E){n{q) I $_,^(0) \n{k)7r{-k) )IT out( 7r(fc )7r(-fc ) | J(0) |0) e'^^"^^)* , (103) 



27T 

where the phase space factor J\f{E) is given by 



Af(B) = |§. (104) 



and the scattering amphtude Ai{2Eq) = M.{q, —q; q', —q') {\q'\ = \q\) is given by eq. (pSf ) 

{2nfS^'^\pm-qi-q2){iM{qi,q2;m)) = ont{m\n{qi)7i{q2))in - ont{m\n{qi)n{q2))ont- 

(105) 

Since J creates a two-pion state of fixed isospin and angular momentum from the vacuum, 
we have 

out(vr(g)7r(-g) | J(0) | 0) = e*^(^^') | out( vr(g )7r(-g ) | J(0) | 0) | . (106) 

On the other hand, the elastic scattering amplitude is given, in terms of the scattering 
phase, by 

M{2E,) = e^^(2^^) sin(5(2E,)) , (107) 

and satisfies the unitarity relation 

2Im[A^(2Eg)] =U{2Eq) \M{2Eq)? . (108) 



Using eqs. ( p.06| )-( pi78D , we obtain 



<3^it,q) = ^e-'^''' |out(7r(g)7r(-g)| J(0)|0)| cos5(2E,) +••• , (109) 



where the ellipses represent the principal value term in eq. (|103| ). 



Appendix C 

In this appendix we reproduce some of the results of ref. in a form suitable for our 
presentation. We recall that we are studying the eigenvalue problem for a two particle 
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system of mass m in the center of mass, subject to a repulsive spherically symmetric 
potential V{r). For simplicity, we assume that V{r) only gives rise to an s-wave scattering 
phase shift where k denotes the center of mass relative momentum of the two particle 

system. We also assume that the potential has a finite range R, outside of which it vanishes 
identically. Our problem is to find the structure of the eigenvalues fc^ and the corresponding 
eigenfunctions of the Schroedinger equation 

(A + k^)^{x) = mV{r)^{x) (110) 

with periodic boundary conditions in a cubic box with sides of length L > i?. In particular 
we are interested in the cubically invariant eigenfunctions with non-zero s-wave projections; 
these are the only ones which can contribute to the sum over intermediate states in eq. 



As a result of the hypothesis on the range of the potential, for R < r < ^, eq. (|110| ) reduces 
to the Helmholtz equation 

(A + A;2)^(x) = . (Ill) 

We start by considering the case ^ ^n'^, where n denotes any vector with integer 
components. In this case there is no solution which is valid throughout the whole volume. 

In order to treat this problem we need the finite volume Green function 

{A + k^)Gv{x - x) = 6{x - x') , (112) 

given by 
where 

p^ = —n. (114) 

It is convenient to consider the expansion of Gv{x) in spherical harmonics. In particular 
we will be interested in the spherical projection of Gv{x) 

G<,°'(r) ^^fdn Gv(x) = i V j^'^ (115) 

iPnf 

with Pn = \pn\- For r 7^ 0, G'y^(r) satisfies the Helmholtz equation, so that we have 

(0) cosfcr sin/cr 
Gy {r) = z h c — . (116) 



'For the treatment of the most general case we refer the reader to ref. 
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The term — in eq. ( |116| ) by itself satisfies: 

(A + fc^)(-^) = 5(f-x') (117) 

so that the second term on the r.h.s. of eq. ( [L16| ) must be regular. A possible way to 
determine c is from the relation 

c = \im-^ (rG'^\r)) (118) 
'•^0 dr 

so that eq. (|116|) gives 

cos 



hm 



;ii9) 



We add a note of caution at this point. The sum ^J2{p„} k^-^ formally the same 
as Gy(0) and is therefore divergent. The limit in eq. (|119|) could provide the necessary 
procedure for the analytic continuation (analogous to s-analytic continuation in ref. 0). 
We will not expand further on this point and we note that eq. (|116|) , together with the 
value of c given in eq. ( |121[ ) below, coincides with eq.(3.29) in the second paper of ref.|^. 



Eq. ( |119| ) together with the definition 



VAtt frf [n^ - r) 

{n} 



gives 

2'oo(l,g^ 
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When R < r < ^, the potential is zero and the wave function propagates according to the 



Helmholtz equation eq. ( |111| ). Then, as usual, knowledge of the Green function Gyix — x') 
implies that 

-^{x') = I d^x 4 ■ [ip{x)d^Gv{x - x') - 4^(f)Gy(f - f )] (122) 

Jv-Sr 

where the integral is extended to the cubic volume with the sphere of radius R subtracted. 
Using Gauss' theorem and the periodicity in V, we have 

^(f ) = R'^ [ da [ipix)d^Gvix- x') - 4V'(^)Gv(x - f )] " n (123) 

JSr 

where n is the internal normal to the sphere Sr. il){x) itself is a very complicated object, the 
superposition of an infinite number contributions from states of different angular momenta. 
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The quantization condition is that its s-wave projection must coincide with the unperturbed 
s-wave wave function in infinite volume: 

„ , , sin fcr ^ cos kr , 
f(r) = A— - B— (124) 

where A and B are related to the s-wave phase shift 5 by 

tan (5= — . (125) 

Here we have used the condition that only the s-wave component is interacting. In general, 
the /-th partial wave contributes a radial function: 

fi{r) = Aai{kr) + Bini{kr) (126) 

where ji{kr) is regular everywhere, while ni{kr) is singular at the origin and signals the 
presence of the interaction 

tan5z = -^. (127) 
Ai 

Since we assume that only the s-wave component is interacting, it is only for the s-wave 



that B 0. Now if we set Bi = on the r.h.s. of eq. (|123|) , the result must be exactly 
zero, otherwise we would find regular solutions of the Helmholtz equation with correct 
boundary conditions, which is impossible if fc^ 7^ for all p„. This means that in this case 
we can put to zero all contributions from states with Z 7^ on the r.h.s. of eq. ( |123[ ). 

Since the only contribution to the r.h.s. of eq. (|1 23|) comes from the rotationally invariant 
component of iIj{x) and we are interested in the projection of eq. ( |123| ) onto the s-wave, we 
require the angular average of the Green function Gy\\x — x'\) with respect to x or x': 

Gp''\r,r') = [ —GP ( Vr^ + r'^-2rr'cos9) . (128) 

Evaluation of the integral gives 

^(0 0)/ /\ sin fclr — r'l — sin fc(r + r') cos fc(r -|- r') — cos A;(r — r') 

G,/ (r, r) = ; c — . (129) 

V K ^ J 87rA;rr' 2k^rr' ^ ' 



The s-wave projection of eq. ( |123|) then gives rise to the quantization condition 

f{r) = -A7:R'[f{R)^G''P'\Ry)-f\R)Gf'\Ry)] with r' > R . (130) 

After some elementary algebraic manipulations eq. (|130|) reduces to 

sin /cr cos kr B( k cos kr — An c sin kr) 
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from which the Liischer quantization condition, eq. (^) follows immediately. 



The derivation sketched here shows that the Liischer quantization formula enumerates all 
the states with k"^ ^ t^"^^ that have a non-zero projection on the s-wave. These states 
are non-degenerate, since the final eigenvalue equation has only one solution. There are 
other eigenstates, not described by eq. (|^), which must be considered; these include those 
with k"^ = for some n. These states have an energy which is independent of the 

interaction because we have assumed that only the s-wave phase shift is non-zero 0. They 
are constructed as follows. We start with a plane wave of the form e^^"'^. This is not a 
solution of eq. ( |110|) because it has an s-wave projection, so that the potential V{r) acts 
on it in non-trivially. However the combination 

xl,{x) = e'P'^-'^ - e'P""^ , (132) 

where p„ 7^ p^, but p„ = p'^ is an exact solution of eq. ( |110| ) with fc^ = ^n^, because 
it has zero projection on the s-wave. Moreover if p^ is different from the cubic rotations 



of Pn, which is possible when \n\ > 8, the state in eq. ( [132| ) can be averaged to get a 
cubically invariant eigenstate. However, it is clear from this derivation that these states 
do not contain an s-wave component and do not therefore contribute to the sum over 



intermediate states in eq. (^6[) . 



There is the further possibility of states corresponding to the eigenvalue k"^ = and 
with a non-zero s-wave projection ||^. The condition for the presence of these states 
can be found using the procedure outlined before. Since the eigenvalue is fixed a priori, 
the corresponding eigenvalue equation represents a consistency condition in the form of a 
relation between n^, the phase shifts and the volume V . These states therefore are not 
generic and can only occur for particular volumes and/or potentials and we therefore do 
not include them in our argument. 

We end this appendix by briefly showing that the knowledge of the finite-volume Green 
function eq. (|113|) and its projection on a given partial wave (such as the s-wave projection 
in eq. (|116|) ) allows one to readily reproduce the class of summation formulae given in 
eq.(2.51) of the second paper of ref. 0. These formulae are central in computing finite- 
volume corrections in perturbation theory in quantum mechanics and field theory. For 
illustration we will discuss the particular case of a spherically symmetrical function /(p^), 
for which we prove below that: 



{pn} 



"'^^In the presence of fully interacting higher partial waves, these states would have a non-trivial quanti- 
zation condition which could be found by the methods of ref. § or by a natural extension of the method 
presented in this appendix. 
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up to exponentially small corrections, provided that f{r), the Fourier tranform of /(p^), 



(2vr) 



oo 



f(r)^-^l rf^pe-^^W), (134) 



oo 



vanishes exponentially with r = In eq. (|133|) pn = \pn\ and fc^ is different from any of 
the allowed pf^ on the given volume, c is defined in eq. ( |121| ) as a function of fc^ and L. To 
prove eq. ( p.33| ) we start with the finite-volume integral: 

d^xf{r)Gvix) = [ d^xf{r)GP{r) 
V Jv 

j3 \ ( cos{kr) sm{kr)\ 

d'^xf{r)( + c ; M . (135) 

\ Attt kr J 

The assumption of the exponential decrease of /(r) at large r allows us to extend the 
integrals in eq. ( |135|) to infinite volume with an exponentially small error in the volume, 
so that: 



J d'xf{r)Gy{x) = ^EfcPlll' (136) 

{Pn} 



Anr J (27r)3 k^ - p"^ 



d'xfirf^^ = m. (138) 



Eq. ( |133| ) can also be used in the limit as k approaches one of the K say. In this limit 



This proves eq. (|133|) . 



both sides of eq. (|133|) become singular and we can write: 



where, as usual, uk is the number of p„'s with pf^ = and the prime on the summation 
indicates the omission of the vk terms with p\ = K^. From eq. ( |121| ) we have that 

where c^'^^ is defined as in eq. ( p.21| ) but with the vk terms with = n\ = {LK/2ttY 
excluded from the sum in eq. ( |120| ). c^^^ is therefore regular in the limit fc^ — ^ K^, 

lim d""^ = ^ , (141) 
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where the number zk is defined by 
Eq. ( |139| ) can now be rewritten: 

{Pn} 

where /'(/c^) represents the derivative of /(A;^) with respect to k"^. Eq. ( |133|) and ( |143|) 
can be easily extended to functions with an angular dependence and to higher order poles. 
We stress also that if we take V 00 keeping the physical momentum K fixed (by 
appropriately tuning the sequence of volumes to make this possible) then we must consider 
Tlx, Zk and vk to be functions of the volume and therefore the finite- volume corrections are 
not simply given by the explicit factors of 1/L and 1/V on the right-hand side of eq. ( |143| ) . 



Appendix D 

In this paper we have seen examples of correlation functions from which one can extract 
the modulus of the amplitude and others from which one can obtain the modulus of the 
amplitude times cos{5{2Eq)). The physical information which can be obtained depends 
on the correlation function which has been computed. In this final appendix we present 
one more example, demonstrating that cos{5{2Eq)) can be determined directly from the 
evaluation of four-point correlation functions. 

We start with the following correlation function 

GA{ti, h, ts, t4) = ( I $,-(ti) $_,-(t2) 4>K^, h) ^\0, U) I 0) , (144) 

with ti > ^2 > ^3 > ^4- As before we assume that we can neglect inelastic contributions 
and that the fiavour quantum numbers of the interpolating operators are chosen so as to 
create two-pion states with a fixed isospin. For large ti and large negative ^4, single-pion 
states dominate in the time intervals (ti,t2) and (^3,^4) so that 

^4(ti,t2,t3,t4) = ^e-^'*^ j ^ ^e-^^l*^l(7r(g)|<|._,^(t2)0(O,t3)|vr(^)) . (145) 

Thus we need to consider the matrix element {7r{q) \ $_g'(t2) 0(0,^3) | vr(fc) ) to which we 
apply a similar procedure to that in appendices A and B. Inserting a complete set of 
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(two-pion) states between the two operators and applying the reduction formulae we find 
that 

^4(ti, h, h, U) = ^ e-^^(*^+*-*3-t4) [1 + LX{2E,)M{2E,)] [l - '-Af{2E,)M*{2E,)] 

+ P.V. 

= -|;^ e-^«(*i+*^-*3-**) cos\5{2Eg)) + P.V. , (146) 

AEg 

where P.V. represents the principle value integral over energies other than 2Eq. Thus from 
the evaluation of the four-pion correlation function and the determination of the coefficient 
of the exponential with exponent proportional to Eg we may obtain the cosine of the phase 
shift directly. 
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